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The inherent symmetry of a spherical reflector microwave antenna 
suggests that if the fields incident upon the reflector from a feed 
required to produce a specified reflected field can be determined, 
then scanning of this reflected field may be accomplished by rotation 
of the feed assembly rather than by rotation of the entire antenna. 

A technique has been determined to solve for the fields in vicinity 

of a hypothetical transmitting feed that will, following reflection 
from the surface of the spherical reflector, yield a specified field 
pattern at the reflector aperture. From given fields over a sphere, 

a portion of which coincides with the surface of the spherical reflec~ 
tor, a spherical harmonic series expansion is used to calculate the 
fields on a smaller concentric sphere encompassing the feed region 

in solution of the boundary value problem. 


This technique is extended to a case where it is shown that 
geometrical optics approximations are valid over a specific region 
of the reflector system, and that by ray tracing, the fields on an 
mapetrary reference spliere, wiose radius is confined only to the 
axis of symmetry of tho reflector, can be determined. From this 
field distribution, the spherical harmonic series expansion can be 
used to determine the fields on the surface of a smaller concentric 
sphere enclosing the hypothetical feed. It is shown that this tech- 
nique has flexibility wherein a number of origins could be examined 
for various feed spheres with fewer terms of the series expansion 
than with the fixed-origin case, 
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1.1 


CHAPTER I 


GENERAL BACKGROUND 


The fact that a reflecting parabeloid ef revelution, when il- 
luminated by a source of spherical waves placed at the geonetrical 
focus of the paraboloid, will preduce a uniform plane wave over the 
eperture of the reflector is well known. The generatien of a pencil~- 
shaped bean antenna pattern has reached a state of refinement such 
that extrenely narrow beamwidth, low sidelobe levels, reasenably 
bread~band operation, high aperture cfficiency and gain are readily 
attainable for many applications. 

The axis of this pencil beam coincides with the axis of revolution 
ef the reflector, and attempts to direct the antenna pattern by off- 
axis rocHont of the point-source Food have failed fer all but the 
smallest offeaxis angles. Therefore, to direct the pattern requires 
that the entire reflector and feed assembly be rotated or moved to ain 
the antenna pattern in ire direction of interest, 

Widespread use of large aperture paraboloidal reflectors with wide- 
angle scanning requirements has stimulated considerable interest in 
developnent ef impreved feed systems to alleviate problems of ground 
noise radiation, spurious pelarization effects end low aperture ef~ 
ficiencies encountered in radio-telescope applications. Other such 
problems that. are encountered include difficulties in maintaining 


mechanical and structural telerances of the reflecor surface at partic~ 


ular frequencies of interest, 








The inherent symmetry of the spherical reflector suggests that it 
could provide some relief from the problems previously described. Assum- 
ing an energy source could be obtained that would, upon reflection from 
the spherical surface, provide an antenna pattern whose beam coincided 
with the radius of the reflector upon which illumination occurred, it 
would be no longer necessary to overcome the inertia of the entire 
system to meet a wide-angle scanning requirement, but merely to direct 
the feed assembly as desired. In such a case, it is feasible that large 
apertures could be illuminated to provide desired gain . Further, the 
symmetry of a spherical section enhances both ease of fabrication and 
verification of system tolerances in many applications. Other advantages 
that may be realized from such a system follow directly from speculation 
of possible application, including for example an antenna system whose 
reflector could be imbedded in the earth, with few moving parts, pre- 
senting a low silhouette to wind forces. 

Unfortunatly, analysis of the spherical reflector reveals its 
inherent spherical aberration. For a plane wave incident upon a 
spherical reflector, rays close to the axis of symmetry intersect at 
a point defined as the paraxial focus, exactly half-way between the 
center of curvature of the spherical section and the vertex. Other rays 
reflected from the surface intersect the axis at points closer to the 
vertex. This phenomenon may be observed in figure 1, where rays are 
constrained by Snell's law to reflect with an angle equal to that of 
incidence, as measured from the surface normal at the point of interest. 


A more complete explanation of this reflection phenomenon will be given, 
7 
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Figure A 


This suggests that rather than a distinct point focus there exists a 
fuzzy focal region. Illumination of the reflector with a point source in 
the focal region reveals a diffraction pattern with wider beamwidth, 
higher sidelobes and much less gain than a paraboloidal reflector of 


identical aperture. It is this aberration that has presented an obstacle 
to development of the sphcrical rather than paraboloidal reflector fer 
use as a microwave antenna. Since intreduction of the problem to the 
microwave community in WorldWar II, however, there has remained an 
active interest in the problem, Mest receatly, with censtructien in 
1960 of the Arecibo (Puerto Rico) Ionospheric Observatory's 1000-foot 
diameter spherical eines considerable interest has bcen generated 
in the solution to problems of spherical aberration at radio frequencies 
by means of various feed systens, 

1.2 It is appropriate to review past efforts relevant to the problen,. Early 
attempts to explain optical phenomena followed from geometrical theory. 


Euclid presented an analysis ef the reflection laws arcund 300 B.C. while 


observing that light energy could be concentrated by means of 
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reflectors. Early astronomers used spherical mirrors as primary 


collectors of light incident upon their telescopes. Presumably the 
ease of manufacture of these spherical mirrors led to their widespread 
use, even though the geometrical analysis of the reflecting propertics 
of paraboloidal surfaces had been conducted. The difficulty of manufacture 
of the latter apparently precluded exploiting thelr point~focus chzrac- 
Geristic. 

The concept of a spherical reflector for microwave use evolved from 
the World War II advances of radar technology. In 1941, Ashmead and 


Pipperd suggested two approaches to the aberration encountered: 


1. For reflectors of large radius and restricted 
aperture, the surface of the spherical reflec- 
tor departs only slichtiy from that of a para 
boloidal reflector. It follows that the diffrac- 
tion pattern is only siightly different from that 
of a paraboloidal reflector. If the slight depar- 
ture from plane wave illumination at the aperture 
and the value of gain accompanying the restric- 
tive aperture can be tolerated, this approach 
permits a simple design and is suitable for wide 


angles of scan, 








Ze For the general case of spherical. reflectors, 
should the primary energy source be a special 
radiator that provides not spherical waves but 
waves whosSe departure from spherical waves is 
just that necessary to compensate for the differ- 
ence between the sphere and paraboloid, then 
reflected wave fronts wil] be plane waves, uniform=- 


ly illuminating the aperture. 


Early work indicated that it was possible to create non-spherical 
waves of approximately the right shape by proper combination of wave~ 
guide horns and dielectric lenses. Unfortunatly, difficulties in 
correcting for the larger phase error did not, at the time, justify 
PuLSswance of this approach. Rather, to meet the expedience of the 
wartime effort, the first approach of toleration of error rather than 
conpensation was used exclusively. Subsequent studies exploiting this 
technique reached a natural conclusion with the work of gts who 
extended and verified the theory of restrictive aperture large radius 
spherical reflectors. His experimental work indicated that with a ten- 
foot aperture hemisphere absolute gain comparable to that of a forty- 
inch paraboloid was obtained over a scan of 140° with rotation of 
the feed at a constant radius from the reflector. 

The second, more interesting approach to the solution allows the 
use of a larger aperture at the expense of increased complexity in 


cost, design and some reduction in the available angle of scan. 
io 











In the context of the latter approach corrective devices have fallen 


into the follewing categorics: 


Ll. HERS structures 
2. auxilary reflectors 
3. line-source feeds 


4, transverse aperture feeds 


Lens structures, while providing satisfactory correction at optical 
frequencies, provide marginal correction in microwave applications, Sizes 
comparable to the reflector and primary source, requirements for large 
pleces of carefully controlled precision-cut dielectric and inherent 
dielectric losses, to say nothing of the cests and weights invelved, 
relegate this category of ccrrective devices at microwave frequencies to 
a pesition of academic intercst. 

Auxilary reflectors are of minor interest when the undesirable charac~ 
teristics ef shadowing er aperture blockage are imminent. Even in cases 
where this is of small concern, the excessive weight ef the secondary 
reflector must be considered, as must the mechanical arrangements for 
support of the assembly. A gregorian reflector has been described by 
Holt and Beuciey wherein a reflecting surface intercepting a great deal 
ef the energy incident from the spherical reflector is installed. The 


Surface of the sub-reflecter is such that it provides a constant path 


length from seme reference plane near the aperture te the primary securce. 
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Experimental results indicate that the scanning capability was realized 
to 30° by rotation of the sub-reflector on a constant radius from the 
spherical surface. Sidelobe levels were very high ee to 
those available with paraboloidal reflectors. While aperture blockage 
was estimated to be 1.5%, no figures on gain weve given. 

The earliest attempts at compensation involved the fact that rays 
incident upon the spherical reflector were reflected to cross the axis 
of symmetry. It was postulated by Spencer, Sletten and eee that the 
incoming energy Peuld be recombined with an appropriate phase shift 
to provide a maximum of received signal. Reciprocity was invoked to 
suggest that a properly phased line source could produce a wavefront 
that would provide uniform illumination over the averture. Techniauecs 
o£ geometrical optics were used for path length calculations. The © 
pemnces considered were open-ended waveguides, polyrod travelling wave 
structures, discrete dipole sources alog the axis and slotted waveguide 
feeds. Indeed, the first feed system constructed at the Arecibo 
reflector was based upon a modification of this Gestion! Aberrationless 
scan of over 60° was predicted for this initial installation~-as well 
as reasonable aperture efficiency. Actual results were well under those 
anticipated. The high-power dual circularly polarized feed illuminated 
the full aperture of 1000 feet with an efficiency of 21% at 430 Miz® 
Love” described a composite radiating system wherein channel guide 
feed and a polyrod array provided a travelling wave structure to illumin- 
ate from both the line focus and the paraxial focus, He predicted a scan 


of 110° for a hemispherical reflector, and efficiency on the order of 
(2 
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55%. Schell Osnvestigated the £ields along the axis by application of 


optics solution in the non-zero wavelength and agreement in the zero~ 
wavelength case. Schell's solution consisted of the use of a sufficient 
number of discrete sources located along the axis, illuminating partic- 
ular regions of the reflector, The analysis was based upon the fields 
present along the focal axis when the reflector was illuminated by an 
incoming plane Wave. In an attempt to solve the problem based upon 
Schell's solution, Cohen and ee constructed a feed system at the 
Arecibo reflector. Thelr compound antenna feed consisted of a four- 
element Jinear array with a reflecting screen and a siztecn-element 
linear array located on the axis. Experimentation determined that gain 


~ 


was lower and sidelobes higher than had been anticipated. Aperture 
efficiencies were on the order of 30-38%, Necornieee followed this 
study with an analysis that considered the transmitting case and, in 
particular, the aspect of polarization. His enalysis concluded that a 
gain deterioration was inevitable with a circular waveguide with cir- 
cumferential slots that did not consider the azimuthal fields as well. 

as the longitudinal fields. The theory of McCormick was used by 

Love and Giuenee to illustrate the feasibility of illuminating from 

a leaky cylindrical waveguide that considered these fields, Mcst 

recent Jiterature indicates that a line source feed with high performance 


has finally been achieved at Arecibo. LaLonde and soceeigee note that 


their feed has been installed at the observatory and has realized an 
13 








overall efficiency of 70Z at 318 MHz. Designing for a flat waveguide line 
feed with transverse slots was necessary to avoid problems of previous 
feeds which, it was felt, did not provide sufficient control over 
amplitude and phase variation along the line focus nor sufficient 
isolation from the effects of mutual coupling between the radiating 
elements of the feed. The lower frequency of operation was specified to 
avoid losses due to surface imperfections of the reflectox. The design 
of the phase and amplitude required was based upon a geometrical optics 
analysis of the rays intersecting the axis. This design is significant 
in that it represents the first installation of a corrective feed in 

a spherical reflector that has yielded aperture efficiencies approaching 
these of pcraboloidal reflectors. Aside from strengthening arguments for 
fixed~-reflector scanning antennes, providing some 40° of sean with high 
SEF iciency at Arecibo, the success of this feed indicates the potential 
of the spherical reflector antenna for fully steerable arrays, where 
desired. 

It has been noted by peeae die thet Line Source radiator correcting 
feed s ens must be capable of producing a wave whose polarization and 
shape axe not easily obtained with elementary radiators.This basic fact 
has caused present day iine source correcting ieee to fall short of 
expectations, with exception cf the LaLonde-Harris feed. It should be 
noted, however, that the requirements of Arecibo call for high power 
with capability for dual circular polarizatid:, The line feed of LaLonde 


and Harris achieved success in providing for linear polarization. While 


the high power requirement may be easily met with the slotted guide, 
14 : 








the authors note that polarization flexibility has yet to be developed. 
Presumably, the transverse aperture corrective feed designs have no 
jmuhepent difficulty in polarization flexibility. Banlyeattempibe to 
provide aperture corrective feeds met with reasonable success in over- 
coming the aberration, but were hampered by aperture blockage of the 
moderately large fced systems. Sletten and Mavi oleae were succesful 

in lowering sidelobes substantially, while plagued with aperture 
blockage. An aperture feed was suggested for the spherical reflector by 
Burrows and Ponda in which they predicted a possible gain of 61 db. 
compared to the maximum gain of such an aperture the size of Arecibo at 
62.8 db. Their feed system was derived from a spherical harmonic expan= 
sion of the fields necessary to produce a uniform distribution over the 
aperture, based upon the fields incident upon a spherical reflector. 
Spencer and ydhae conducted a eeeiicd study of the focal region ficlds 
associated with the spherkeal reflector illuminated by incoming plane 
waves and concluded that only a vector field solution considering polar= 
ization would yield adequate information, They solved the diffraction 
intesral by means of a stationary phase approximation, noting the 
Similarities to the geometrical optics approximations, Their theory 

was shown to have reasonable agreement with fields actually found in the 
focal region by experimentation, and their work has been used as a 

basis for further study into aperture feat: Their work, however, did 
not offer a solution to the problem, but rather provided a theoretical 
analysis of the fields that could be expected. Work by Minnet and 


, 22 : : 
MacA,. Thomas has been extended to show that the fields in the image 
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of synmetrical focussing reflectors can be generated by certain hybrid 

“ AY Gates: : Ad m4 2 3 r 
modes to provide theoretically aberrationless behavior ~. An aperture 
source iswdescribed that uses a corrugated wavegusmder"strictureto generate 
a set of waves that combines to provide a theoretically high efficiency. 
Further developing his earlier ideas of spherical harmonic expansions, 
Ricardi has recently provided a technique for synthesis of the ficlds 

ee a4 

on @ surface enclosing a transmitting feed. This treatment for the 
transmitting feed synthesis provides for the solution to an electro- 
magnetic boundary value problem, identifying as boundaries the reflector 
surface and the surface of a smaller sphere enclosing the feed region. 
It is with this technique, substantially different fron those involving 
geometrical optics or approximations of stationary phase, offering a 
solution that appears to contain inherent accuracy, with which this 


thesis is concerned. 
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CHAPTER IT 


PROBLEM STATEMENT AND THE THEORY OI GEOMETRICAL OPTICS 


2.1 The solution propesed by Ricardi te the preblen of the spherical 
reflecter microwave antcuna relatos twe sets of clectromagnetic fields 
ever portions of surfaces of two concentric spheres. Due te the mathe- 
matical formulation of the problem, one surface coincides with that of 
the reflector while the ether surface is taken at a radius of interest 
Within that ef the reflector. Origin of both spheres is at the center 
ef curvature of the reflecter (figure 2). The fields over the reflector 
surface are specified te be those fields that will produce a desired 
frela distributzen ever the aperture of the reflector felleving reflection 
from its surface. The fields at the reflector surface are then exnanded 
in a spherical wave expansion. Thus it may be determined that if the field 
distributicn given by this expansicn, evaluated at the inner surface, can 
be somehow generated by a source within this surface that the appropriate 
reflected fields will appear, as desired, over the aperture, Interpreta- 
tion of the fields at various radii has revealed that there are certain 
radii coincident with nearly constant phase fronts. Synthesis of these 
fields would be enhanced due to a potential lessening of the conplexity 
ef the radiating structure. Exactly which radius to be chosen for the 
inner surface is beyond the scope ef this work, but is covered by 


Pileardi-- 








It is the hypothesis of this work that the technique of solution 
to the boundary-value problem may be applied to a reference sphere 
whose origin is a variable distance along the axis of symmetry and to 
the surface of which outgoing rays from ea hypothetical feed located 
within this reference sphere would be traced if the spherical reflector 
were absent. The rays outgoing from the region of the hypothetical feed 
would be precisely those rays representing the fields that would, upon 
reaction with the surface of the spherical reflector, reflect a desired 
ililumination over the reflector aperture. These rays would be based upon 
the rays reflecting from the spherical surface when illuminated by an 
incident plane wave over the aperture. Strictly speaking, this analysis 
begins with the principles of the receiving case, but then requires 
the field distribution to be produced over the reference sphere ~ some 
transmitting feed. It is then the reaction of the radiated field with 


2 


the spherical reflector that theoretically produces a desired wavefront 
over the aperture, The geometry for this case is given by figure 3. 
Thus, within the limits imposed on ray tracing technique by theory 
of geometrical optics, Ricardi's results for the transmitting feed 
case may be compared to an approximate solution that originates with 
Mae recciving situation. It is further possible that Such a solution, 
known to be strictly valid only in the zero-wavelength limit, may yield 
results with less computation time than the exact solution within an 
appreciably small deviation in results. The choice of origin of the 


system of concentric spheres sugeests that a more detailed analysis of 
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the fields in a particular region may be conducted in an ccononical 
fashion, since the fields will be distributed over a sphere of smaller 
radius than the previous problem's inner sphere with origin at that 

ee the reflector. 

We thus proceed to describe the illumination of a spherical. 
reflector by a uniform plane wave to illustrate certain character- 
istics of the reflector. This is followed by development of geometrical 
optics principles eseancial to the deseription of the fields over 
the reference sphere. Finally, expressions will be developed that 


yield the amplitude, phase, and direction of the electric field 


over the reference sphere. 


2.2 It has been noted previously that rays reflected from the spher- 
ical surface following illumination by a uniform plane wave over the 
aperture cross the axis of symmetry at varying distances from the 
memaxial focus. This phenomenon may be analytically formulated to 
illustrate the generation of caustic surfaces by these reflected rays. 
In fact, it will be shown that the caustic surface is the envelope 


of reflected rays. Consider the situation given by figure 4. 











Figure G 


cquslic viewed (n plane oe Tees 


Constrained only by the fact that the angles of reflectien and incidence 


are identical, the equaticns of rays may be given by the slope intercept 


form of the straight line equation. That the reys are straicht lines and 
2 ro) 


the angles equal will be covered in detail in latcr sections, Viewing 


a 


the system in the x=-z plane we have: 





Points on the reflector surface are given by: 


x = R sind 


z= R cos8 


(Zat) 


Letting points on the envejonve he xee Zs and substituting into 2,1: 


fa 








Xo = R sind - (R cosd ~ zd tan20 Caz 


Thus points on the caustic surface as viewed in the x-z plane become 
functions of parameter ©. We then form the function f£(x,z;0) = O that 


will be satisfied by points on the caustic. 


E(xy zo) = ee R sin® + (R cosd = z )tan20 ddan Sy) 


Adjacent points on the caustic generated by rays incident at angle 

eeoo Will still satisfy the function for small 60. Since both functions 
£(x 5%) and E(x 9% 0160) are equal to zero, their difference as well 
mequal to zero. We can then form the expression valid in the limit of 


Small 60: 


. 











E(x 42 0160) = £( 42,9) _ 
liméos0 ——S—S cs 6 0 
50 
or def (% 452 90) 
Gc 
“ Se d 
46 0 (2.4) 
Differentiating and solving for 2 viieldss 
cr eae ae sin20 + 1)cosé G5) 
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Zc. 








This expression for z is then substituted into 2.2 to solve for xO 
Manipulation with trigonometric identities yields the parametric form 


of the equations for the caustic surface: 


sce R sind = R( 3 sinO ~ sin30) (226) 
d, 


Z.°5 Reos ( 2 sin?9 + 1) = R( 3 cosO ~ cos39) (2.7) 


R 
2 4, 
Thus it has been shown that the reflected rays‘ envelope forms the 
caustic surface and quantitative expressions for this surface have been 
developed, They will be called upon in conjunction with development of 


the expression for amplitude. 


» 


pea esectromagnetic field problems may be formulated in two contrasting 
fashions, following from a choice of application of the Maxwell field 
equations. In one treatment, relationships may be derived that determine 
what fields will arise from a prescribed set of sources. The field and 
vector Helmholtz equations are integrated by application of a vector 
Green's theorem, and express the fields at an observation point as the 
sum of contributions from the sources distributed throughout a particular 
volume and from fields on the surface of the volume arising from sources 
outside the volume.” This rigorous treatment indicates that the field 


represents a flow of energy outward from the region of the sources, and 


25 








that these sources and fields satisfy certain assumptions of continuity 
over the surface of interest, as required by application of the Green's 
theorem. In the second treatment, the integral relations derived by thea 
first method may be applied to yield expressions for the field vectors 

at a specified point,egiven the values of the fields over a surrounding 
equiphase surface, or wavefront, without direct reference to the sources 
generating the wavefronta/ These integral expressions are the analytic 
formulation of the Huygens-Fresnel principle, which states that each point 
on the ecquiphase surface can be regarded as a secondary source of elec- 
tric and magnetic current and charge which in turn gives rise to a spher- 
ical wavelet. The fields at a point of interest may be obtained by a 
superposition of these wavelcts, considering the phase difference upon 
arrival at the point of interest. Hence, energy flow is a wave phenomenon- 
ane which may be studied without reference to the sources of cnergy. As 
before, the rigor of the solution depends upon adherence of the given 
field distribution to several constraints as imposed by application of the 
Green's theoren. 

While the Huygsens-Fresnel relation provides expressions for solutions 
to the wave equation that satisfy Maxwell's equations, complexities of 
integration of these expressions often render solutions impractical, or 
require approximations to be made that eventually lessen the rigor of the 
solution only after substantial manipulation. Jn these cases it becomes 


convenient to approach the propagation of waves from the standpoint of 


geometrical optics, wherein a rigorous solution is sacrificed for a 


AG 








simpler formulation that offers a more immediate solution, subject to 
several additional constraints on the problem. In geometrical optics, 
successive positions of cauiphase wavefronts of the fields and an 


associated system of rays are related. 


In an arbitrary medium the wave field is characterized by both a ray 
velocity and wave velocity at every point, where the ray velocity is the 
velocity of energy propagation, directed as the ray passing through the 
point of interest, and the wave velocity is the rate of displacement of 
the wavefront in the direction normal to the surface. The energy propa- 
gates at the speed of Jight in vacuum, c. The wavefront propagates at the 
velocity v, a function of the mediun. = 

We proceed to describe the relations among wavefront and phase, con- 
sidering wavefronts generated by some source propagating in a steady- 
state radiation. Let the wavefront at time to be the surface LOsy,Z)*L,, 
a surface of constant phase ¥ = (w/e) L relative to some point. The wave~ 
front at some short time 6t lates is the surface L(x,y,z2) = Lo + OLs 
The phase difference between successive wavefronts is then given by 
(w/c)6L. Since the wave proceeds from one surface to the next in time ét 
while the phase at any fixed position changes at the rate w, this phase 
difference must be equal to wét. Further, if oer is the distance between 


the adjacent surfaces and if v is the wave velocity, then vét = oe: 


25 








We then relate these expressions: 


(w/c)é6L = wét =(u6 _d/v 


or 61/6. = c/v =n , index of refraction of the medium (9.8) 


Note that expression 2.8 is the directional derivative of a scalar L 
ap the direction of the normal to the wavefront. By definition of the 


directional derivative of a scalar function ioe 


dé/ds = T+V6 


where T is the unit vector indicating direction in which the derivative 


is taken, By definition of the vector dot product, with 0 the angle 


> 


between T and V9, 
do/ds = TeV¢ = [Tt] | Vo] cos® = |Vo|cosd 


For our case the direction of interest in normal to the front, and O= 0. 


Hence; 


c/v= n (2.9) 


tt 


Ok ee 
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Considering the curvature of rays in an inhomogeneous medium we can 


demonstrate that the rays ere rectilinear in a homogencous medium. Let 


¢ 


Seecta unit vector in the direction of ‘a. nay at the pate or apes . 


This vector is normal to the wavefront and has the directien of VL. Then, 
MOE) Sa ere alive | im = 1 
Then Se> Vina (2.10) 


Let N be e unit vector in the direction of the radius of curveture ef the 


Yay at the same point of interest and p the radius of curvature of the 


29 


ray at the point. The vector curvature is defined” as N/p . This vector 


curvature is also given by dS/ds, where s is a distance measured aleneg 
the ray. By definition and vector identity the following expansicn can 


CY 


be made:3 


i 


dS/ds = (S*V)S Horas ome 11/0 (2,11) 


This expression can be manipulated to show a simple relaticnshin between 
curvature of the rays and index of refraction of the inedium. Details are 


left to appendix J]. The result of interest is: 


1/p = NeV(in n) (2512) 
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Thus for a homogeneous medium where the index of refraction is constant, 
or independent of position, we see that the right hand side of 1.12 

is zero, indicating that the radius of curvature of the ray is infinite, 
indicating that the rays are,, in fact, rectilinear. It follows that 

VxXS = 0, which is a sufficient condition for the existence of a family 


of surfaces orthogonal to a field of vectors S. 


The relationship of the amplitude of fields on successive wavefronts 
can be determined from the premise of geometrical optica that rays are 
lines of flow energy and that no power flows across the sides of a tube 

e @ ° q OQ. 4 30 
of rays. These results were verified in the preceeding paragraph. Silver 
considers two surfaces of waves L, and L, and a tube of rays that cuts 


out clements of area dA and dA, on the respective surfaces. By con- 
scrvation of energy requirements, the flow across any section normal to 


the tube will be constant. In terms of the Poynting vector defined by 


a 


* 


S = i(EXH ) 


we can formulate the flow constraint by 


since in free space the Poynting vector is proportional to the magnitude 
squared of the ficld, between successive wavefronts the following relation 


@2ii1 hold: 


Zo 








Eee Ente ace (2.13) 


Considering the ray through point A on the surface Ly coincident with 


the z-axis in figure 5, let the xz,yz planes coincide with the principle 


planes of Ly at A, A ray through an adjacent point B lying in Ly and the 


wee plane will intersect the ray through A at the point 0.» a distance 


Ry away. This distance Ry thus is one of the principle radii of Ly at A, 


A similar analysis reveals the definition of the other radii Ry based 


upon intersection of rays passing through adjacent points in the x-y 


t 


mine, Lhe point A’ on L, ilies on the ray through point A. At A” Jocal 


Z 


t 


coordinates x", y' are constructed. Since the fronts are spaced apart 


by the distance P, the two principle radii of curvature of the wavefront 


L, at the point A‘ are Ry ~’ P and Ry + P, 


Considering the elements of area dA, and dA. around the points A, A’ 


bound by the curve C, C' on surfaces Ly and L, respectively, we can 


write: 


dA, = f(xdy-ydx) 


Z 
CZ ta) 
dA, = LCxidy'=y'dx") 


From inspection of figure 5, and by the law of sines we can relate 


wWersious distances: 


(a 








aw 
-\ 





Ro aR (715) 


Prerelating dx, dx’, dy. dy’ and substituting into 1.14 we obtain a 


relationship between the area clements: 


dA, = dA, (R, + P)(R. + P) (2.16) 
2 es re 
RR, 


pupSstitution of 9.l6éeinto 9.13 we obtain: 


3, |*dA, = 2dA. (R, + P)(R, + P 
1 2 
For the case of a wavefront diverging from a region of sources this 
expression relates the amplitude at the wavefront farther from the 


source: 
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(2.17) 
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From the preceeding development it would appear that geometrical 
optics is based on the local behavior of the wavefront as a plane 
wave. Indeed, this is the basic tenet of the theory. The plane wave 
solution to the wave equation may be used to show that in the limit of 
decreasing wavelength the geometrical optics solution satisfies the 
Maxwell fields equations. The vector solutions appropriate to describe 
Renee iy polarized clectromagnetic ficlds that behave locally like 


plane waves have the form: 


ma) ING esene oi (otek OL (xsy52)) 
‘ (Zee) 
Fre, 
H = B(X,¥,5Z) oJ (ut kK LCs ¥52)) 


The amplitude vectors may be complex, but phases must be independent of 
position to fit the plane wave solutions, Substituting these forms of 


solution into the homogeneous forms of Maxwell's equations given by 


VXE+ jouH = 0 
C219) 
—juck +V X H= 0 
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and expanding the resulting expressions, these may be solved for the 


vector amplitudes A and B : 


A = oie Vile" DB) aC Vee 
WE jue 


(2.205 
B = ko Vis A) = Ce) 


or, 


Wy jour 


where Ko is the free space propagation constant, Ke = Ou a. ; 
Equations 2.20 may be solved to yicld an expression in terms of A alone. 


Using the definition of ko and the index of refraction 
3 
n= c/v = (ep/e uo) 


we obtain the following expression: 


= 


peepee | VL x VL x A) 1 Cvs Vv x A+ VY x Vix AD + 
ie Tee 
O 


oC Vx V xe) 


ee. C252) 
n ke 


Expanding 2.21 with the vector. identity 
VL x VL x A = LOA) —eaGee VL) 


we ovtain for A : 
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A= sl { VL@AeVL) = AlvL|7}> > 1 Grin v= eee) 
ce jk _n? 
O 


same ey ocr (7E22) 


A similar expression results for B, replacing A with B in 2.22. If 

L and the first and second derivatives of A and B are finite, the last 
Baomcerms of 2.22 are of the order 1/k and 1/k compared to the first 
term. As the wavelensth approaches zero, SA approaches infinity, and 


we are left with the following expressions: 


As =) {VL (AcvL) ~ Alvi? } 
n¢ 
(2723) 


B= ~l_ {VL(BeVL) - B/VL|? } 
n¢ 


For these expressions to yield identity of amplitude of the vectors 
A and B we see that the following must be true for the zero wavelength 


case} 
AeVL = 0 
BeVL = 0 


These two conditions note that A and B are transverse to the gradient of 
L, or that they both lie in a plane transverse to the direction of propa- 


gation. Furthermore, since 2.20 expresses B in terms of A we have; 


Be, 








B= Vc/i{(VL x A) -_1VxA } Crom 


ee SS 


n jkon 


We sat that in the zero vavelength case the second term vanishes 
and we are left awe, én expression that implics vector amplitude B is 
perpendicular to A as well as transverse to the direction of propagation, 
since VL/n is a unit vector in the direction of propagation. 

Thus the Asi vectors of geometrical optics possess the properties 
of plane waves in what might be called the far-zone fields, indicating 
preat distances from the source with respect to wavelength. The condition 
for neglecting terms of the order 1/k, and 1/k for short wavelengths 
require that all associated derivatives are finite. In the neighborhood 
Of geometrical focal points or caustic surfaces, the function L varies 
rapidly, and geomctrical optics predicts infinite amplitudes , which 
must be rejected on non=-physical grounds, At the region of seometrical 
shadow the approximations begin to predict erroneous results due to 
the rapid variation of amplitude over distances small eerie to 
wavelength. In these regions the simple methods do not work, and dif- 
fraction end scattering thcory must be used. 

Further work Et has determined that the Huygens-Fresnel relation in 
the zero-wavelength case approaches the point-to-point amplitude relation= 
ship predicted by the geometrical optics concept of energy flow in tubes 


of rays. Agreement holds in the far-zone approximation, at distances 


greater then wavelength, and reveals the constraint on amplitude 
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A's normal derivative that, besides being finite everywhere, 


LBA << 20 


ow: 


A an x 


This is satisfied in the limit of zero wavelength providing the frac- 
tional change in amplitude over distance equal to that of a wavelength 


is smaJl1 compared with unity. 


Since the basic tenet of geometrical optics requires that the wave- 
front behave locally like an infinite plane wave, scattering phenoinenon 
may be interpreted geometrically, and appropriate laws of reflection 
derived from the ficld equations and boundary conditions. That these 
results are valid can be developed from other principles of sie” 
including Fermat's principle and the Jaw of the optical path, demonstrat- 
ing the properties of optical rays and the reflection phenomenon, or 
Snell's laws. 

We consider the two-dimensional analysis of an infinite plane 
wave incident upon an infinitely conducting ( that is to say perfectly 


reflecting) sheet. The geometry of the problem is described by figure 6. 
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The Atncident fields will have the form: 


ee a ikoz coso. }- jk ox Sind, 
i y i 


ae 


>» 


(2.255) 


~jkoz cosa, + jk ox sind, 


ea 
tt 


c(i, h, cos, +i h sind, ) e 


z 


fi 1 
where ~ is the wave impedance of free space, equal to Gey 


The reflected fields will have the following form; 


| id S ale i 
Te ee jk oz co oe jk x sind. 
r y 
(2.26) 


t- sQ + 4k x si 
E = pied h_cosO. +i h sinO )e ays ie a ore Ws ey 
xy Fs r Zt Y 
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Since at z= 0 for all x and y the boundary condition on the tangential 


ey ; a I “quire: é ee aiaeele = 
eomponent of electric tied requires that (E. ue eae 0 


Then, equating reflected and incident fields appropriately we obtain: 


ec. =e 
i rs 


Thus we note that the following boundary conditions exist at the inter- 


face: 


(Es i Eee os ; 
C2?) 
CE) or ~ ae 


These expressions may be given in terms of a vector relationship: 


(2.28) 


We can examine this vector relationship to determine an expression for 
the field reflected from a Surface in terms of the incident field and 


the surface normal- 4ssuming that the field acts locally as a plane wave, 
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With the vector form of the boundary conditions wa form the cross product 
of the normal n and the first expression of 2.28. We then multiply the 
second of 2.28 with the normal vector n and add the two expressions. 


This yields: 


nxE xn-k (neE )n = -n x E, x n + (nek,)n 
r re af i 


Expanding by the vector identity 


AxBxC = (A*C)B = (A*B)C 


we ocotain the expression 


E = ~E, + 2(n°E,)n (2.29) 
This may be interpreted to express quantitatively the fact that the 
reflected ficld experiences reversal of its sense of tangential 
components in order to provide a zero tangential field at the surface 


er the reflector. 


Thus we have demonstrated that the approximate methods of geometrical 
optics involving point-to-point transformations along rectilinear rays 


between successive wavefronts or surfaces of constant phase can provide 
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adequate phase and amplitude information providing certain constraints 
are observed. Furthermore, thet these fields satisfy Maxwell's fields 
equations and the wave equation within the limits of certain constraints 


has been demonstratec by a zero-wavelensth limit. We summarise the results 


ef e@ometrical optics in which we aré interested as relevant to thas werk: 


1. Energy passing through consecutive cross= 
sections of a tube of neighboring rays is 


constant. 
2. Rays in a homogeneous medium ere rectilinear. 


3. At a surface of vnefleéction, an incident 7 
wavefront behaves locally like a plane 
wave, and the reflection of rays may be 
geometrically interpreted in accordance 


with specular reflection laws. 


The following constraints must be observed for the geometrical optics 


approximations to hold: 


I. All lengths of interest, radii of curvature, 
for example, must be large compared to a 


wavelength. 
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3. 


Spacing between neighboring rays must 
change very little over distances on the 


order of wavelength. 


The fractional change in amplitude of 
field quantities over distances on the 
order of wavelength must be small com- 


pared to unity. 


Geometrical optics cannot be used in 


the region of geometrical shadow, 
Geometrical optics fails at focal points 


or caustic surfaces, predicting non- 


physical results. 


ho 





fee it remains to utilize the previously develop@d expressions to Gesempe 
the phase, amplitude and direction of the electric field over the ref~ 


erence sphere. The choice of radius of the reference sphere is discussed, 


as is the technique of ray tracing to be used. Expressions for phase 


~= 
Nat 


and amplitude on the reference sphere are developed, relating these 
quantities to those incident upon the reflector surface. Finally, the 
distribution over the reference sphere is related to a coordinate systen 
with coincident origin. 

For the geometry of figure 7, we consider determination of the fields 
produced by some source, directed towards the reflector,that will produce 
illumination over the aperture in the form of waves propagating in the 
~ z direction. The fields propagating in this direction will be of the 


form: 


or in spherical coordinates: 


E, H « oJKR cosQ 


Since we desire to specify these fields over a reference sphere of radius 
R' with oricin at z = cR, whére c is some constant fixing the origin, 


the phase expression must be corrected to account for the additional 


vi 
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math Lenoth of travel to the reference sphere, Inspection oo. figure 7 
indicates that the phase variation has two cases of interest, one for 


angi.cs 9 > 0. aid the otherpior theseasaacr 


@D 


< where is the 
oF 5 oe 
anele at which the reference sphere intersects the spherical reflector. 
The variation in field expressions is introduced as? 
kk (Reogaiea 1 
Eu < a oO ( a > 
(27.30) 


( S$ « 
8 < 8 eg tl Sac oe sisele Lo) 


Inspection of the exponent suggests that Lf appropriate values of other 
paraneters could be found so as to keep the difference to ea minimum that 
synthesis might be enhanced, since the smaller exponent Semen 2 
slower varying sinusoid. Furthermore, amplitude considerations in the 
geometrical optics approximation, nanely the desire to avoid the region 
of the caustic surface, suggest that a wise choice would be to consider 
the reference sphere intersecting the reflector surface at the latter's 
extreme edge. Thus the radius of the reference sphere will be constrained 
by the radius of the reflector under investigation, R, the position of 
the origin of the reference sphere’ z= cR, and the maximum angle sub- 


tended by the reflector. By laws of trigonometry we observe: 


12 2 pee one oe 
R Re ac-kK 2Chece ae (2.31) 
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The vay tracing technique procecds from a ray incident upon the 
reflector at a particular angle of interest. The ray reflects in ac- 
cordance with the previously developed lavs of reflection. We then 
assume that if a hypothetical source could produce the appropriate 
fields that would give the same illumination as we have upon the aper- 
ture, these rays would remain as we have determined, but the sense of 
propagation of energy along the rays would be reversed. We then remove 
the reflector, and trace the rays until they intersect the reference 
sphere, Within the Temes of geometrical optics, we can estimate the 
amplitude and phase at each point of intersection, based upon the 
extra path length travelled by the waves, along the rays, to the refer- 
ence sphere. 

It is necessary to calculate the distance L. We proceed with on 
Paetion determining the distance between points on the reflector and 
points on the reference sphere, as measured along the ray of interest. 


Referring to figure 8; 


= ne ats ie es age y) 

My (x. x) on Y > r (2, z) (2732) 
where 

x + R* sind’ cosd' x = R sind cosd 

a m 

es St a te F _ i 

Yan Kk" sin8 ind see R sind sind 

z = R' cosd' + cR z = R cosé 

a m 








Noting from symmetry thet 6 = 6", and substituting, expanding terms 


we obtein: 


ES aE ET 
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2.33) ae (R'sinod' - Rsind)* + (R'cos0' + cR =~ Reosé )4 


To procead, we must relate 0 and 0° . Viewing in the x-z plane, the 


equation of ray is given by: 


x = (z - Zz) tan26 


(2.30) 
0 


) 
Ki 


From the equation of a spherical surface ;: 


R'2 = x2 + (2 =~ ecR)Z + y2 (2,35) 


We can solve 2.34 and 2.35 in the plane y = 0 to yield an expression 
representing the intersection of a ray and the surface of the reference 


sphere. Solving for Zoe intercept of the z-axis we obtain: 


R os Zz 
sin(1- 26) sind (27235) 


Substituting 2.34 into 2.35, and expanding in a straightforward fashion 


we obtain, with substitution of 2.36: 


NS 








g2 =» 22( RsinOsinZd + cReos*26 ) + R2sin26 - cos#20(R'? = c2R2) = 0 


Details of the menipulation are included in appendix J]. Normalizing all 


distances to R, radius of the reflector under investigation, we obtain: 


Ze Zz siné@sin26 + ¢ cos#26 ) + sin20- cos*20(R'? ~ c*) = Q 
OSI), 


where 


R R 


Solving the quadratic 2.37 by usual methods, we obtain for z3 


Zz = ( sinOsin26 +c cos*26) ‘- 


ee ee St 


(sinOsin260 ++ ¢ cos*%20) - sin76+ cos*20(R'? -~ ¢*) 





=< 








(2.38) 
Noting from the geometry of figure 8: 
R'cos ' = 2' & x = eR 
or 
cos6’ = (2 = cR) = (Zz. - ¢c) (2.39) 
—e Fe 
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Substituting equation 2.38 into 2.39 the following relaticnship between 


angles 6 and 6' is cbtained: 


cose’ = .]_ ( sinOsin26 + ¢ cos*20 = c¢ 
“a 


1 








PS 


+ '(sinésin26 + c cos?26 ) = sin*6 + cos*20(R'* ice) ) 








eee 


(2500) 


Equation 2.40 states that there is a relationship between angles such 
that for every angle of incidence on the reflector surface there cxists 
an angle 6', as measured from the origin of the reference sphere, that 
describes the intersection of the reflected ray traced back to the refer 
=e Sphere and the surface of this sphere. Thus the angles may be 
related in terms of the parameters of the reflector under investigation, 
R and 0. » and the center of the reference sphere implicd by the value 
of c. Once the angles associated with a particular incident ray have been 
related, they may be substituted into 2.33 to yield the distance L as 
a function of the point of interest on the reference sphere, 

Equation 2.17 was developed to express the amplitude of field quanti- 
ties on one wavefront with respect to the field quantities on another 
front. It is seen that the amplitude varies with the two principle radil 


of the wavefront and the distance of separation between fronts. Jt is 


necessary, then, to develop expressions for the radii. 
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Once again we consider the case of a spherical reflector iiluminated 
by an incoming plane wave. A typical reflected wavefront has been 
constructed in fipure 9, where the front represents a surface of 
constant phase as defined by the total path length measured along 
@€ach ray, Ly af Jay = R .,. The wavefront is formed by connecting normals 
erected on each rey at the point where total path length equals the 
radius of the reflector. Figure 9 shows the front as viewed in the 
plane y = 0. From symmetry, it may be seen that the surface of the 
front can be generated by revolution of the projection in figure 9 
around the axis. 

At any point of interest, a wavefront may be passed through this 
point by choice of the appropriate path length, At @ach point the 
wavefront will have two principle radii of curvature, Since a ray is 
normal to the front and any two adjacent rays Intersect at the caustic 
surface, forming the envelope that is the caustic surface, it is seen 
that one redius of interest is the distance from the front at the point 
of interest to the point on the caustic, as measured along the ray. The 
second radius of curvature may be determined from the axial symmetry of 
the system. In the $ direction, adjacent rays will intersect at the:-axis 
of symmetry, and the second rere of curvature will be the distance 
from the front at the point of interest to the axis, measured along the 


ray. 
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Rays incident upen the refleeeeh et an angled Ree yellectece. Bits 
anple, pass tangent te the caustic surface and cress the axis, Tne 
ccordinates of the point to which tangency occurs are given by cquations 


2.6 and 2.7, the parametric @quations of the caustic. Fron figure 9 we 





Secs 
=e Sa ee enter eee ne 
Ry MU V(z Z) (x xO 2.4) 
Sey BE: (2, ~ z)° - cam - 0) OD) 
where Zz = Rceosé 
m 
x = Rsind 
mM ne 
Pig rae Sind 
a caeenes 
sin26 


Substituting into equations 2.41 and 2.42, and normalizing all quantities 


to the radius R we obtain: 








Ry =e (2.43) 
2 
7 1 
Ron - omar (et) 
2 cosé 


Detaiis of the manipulaticn involved are covered in appendix J]. 


So 








Although the wavefront has been avalyzed by considcring iliumination of 
the reflector by an incoming wave, it is the fact that if a source can 
be provided that can generate such a wavefront as has been described, tk 
wavefront will be reflected from the spherical surface to yield an 
outgoing plane wave over the reflector aperture. Figure 9 represents a 
graphica], illustration of the method of compensation sugeested by 
Ashnead and Peaparele ee wavefront whose departure from sphericity is 
just that needed to transform into a plane wave upon reflection. 

Since we wish to describe the ficlds over the reference sphere as 
a function of coordinates with origin at that of the refcrence sphere, 
we must introduce a transformation of the field vectors. Initially, 


field quantities were described upon incidence in terms of coordinates 


with origin at that of the reflector. The geometry of interest 1s 
described in figures 10 and 11, considering both cases where the center 
of curvature of the reference sphere is exterior and interior to the 


intersection of the reflected marginal ray (edge of reflector being the 


point of incidence) and the axis. Note that due to axial symmetry there 


is no modification necessary to the 6 components By inspection: 


1, = i. cos (6' = 0) +i. sin (6' - 6) 


r 0 


(2545) 


a Wes : Ee dat See 
toe i. sin (0 6) + i, cos (6 0) 
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Solving for the unprimed components in terms of the primed systen: 


i #i ,cos(! -"6) ~ i 


7 i psin(o’ = 6) 


0 
(2.46) 


i, = i sein(0' - @) + i) ,cos(6' ~ 0) 


From the prececding derivations it may be seen that equations 

tay, 2-30 and 2.29 provide sufficient information to form an expression 
for the geometrical optics amplitude, phase and direction of the electric 
field distribution over the reference sphere, based upon a desired 
reflected field E.. « Quantities of interest, including the radii of curv- 
ature of the wavefront, distance from tne point of the reference sphere 
under consideration to the corresponding point on the reflector, and the 
relationship between angles are given by 2.43, 2.44, 2.33 and 2.40 
respectively. Finally, transformation of the vector components given by 
2.46 relates this field distribution to coordinates at the center of the 
reference sphere. Thus the geometrical optics field is described over the 


reference sphere as a function of the reflector surface under investiga» 


tion and the particuler origin of interest of the reference sphere, 
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In summary, we formulate an expression for the general case where 


a desired electric vector field is specified over the aperture, Choosing 


this field to have wiity amplitude, we describe the field that would 


occur over the refereuce sphere; 


Amplitude: E_ , = F Sion 
rererence het lector | ———= eee 
sphere (R, p ea Be L) 


dk RC cosé + L) 


Phase: Phase( E ) 
reference 
sphere 
Direction: E P z3 ~ - + vee n eE yn 
reference re r 
sphere 
Combining: 
R, RB 4 jk R(cos - L ; 
—— % a: ‘In ‘In - as oO ( ~) { 
refereice enero m reenter ces 
sphere R + ZL R a), 
- ( dn Os 2n ae 


3(2.17) 


(2.30) 


(2.29) 


ook 2(n°E)n} 


(2.47) 








Where: 








Ry = H £988. (2.643) 
So = J) (2.65)) 
Z2cosé 
RO =  j + ¢4 + 2e cose, Os 
L = Te. - R! sind?) + (RE ert + C = COSO)a(2.54) 
cose’ = 1 ( ginOsin26 + c cos*20 ~-c¢ 


an or See 


R! 
n 


“|. emer + ¢ cos?26)%~ sin*6+ cos*26 (Ri? - ¢*) ) (2.40) 


i. = i 1cos (6! - 6) - i.,sin(e' - 6) 


0 
(2.46) 


ee i sin’! - 6) + i,1cos (6! - 0) 
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B42 In this section, we apply the results of the general expression 2.47 
to & particular case of interest. Specifying the field at the aperture 
shown in figure 12 to be a unifoim plane wave linearly polarized in the 


—_ —- 


x~- direction we have: 


se = ee 
aperture x 


Expanding in spherical coordinates; 


; jk_Reosd 
= e oO 


_ L i a 
aperture sin¢ ) (2.86) 


( i .sint cos? + i,cos9cosd - i 


0 ? 


By applying equation 2.29, reversing the sense of tangential components, 
we can determine the sense of direction of the field over the reference 


sphere: 


S6 








direction( FE = i sinOcosd ~ i.cos0cosd -- i. sind REELS, 
Cee eerence r ae 9 : o ¢ ) 


sphere 


Applying equations 2.46 to 2.49 the direction of the ficld vector 


can be expressed in terms of the origin at that of the reference sphere: 


i. zz i 1cos (0! -~ 0) = i,1sin(6! - 6) 


(2 6 46) 
i r: i 1sin(6' - 6) + i,1cos (0! - 0) 


Substituting: 


direction( E ) 
reference 


1 1cosd ( sin€cos(6' - 6) - cos@sin(6' = 6)) 
sphere 


= 


- i,.cosd( sinOésin(0' — 0) + cos@cos(0@' ~ 6)) 


v 


smd beers e hy Yo, 


4 (2.50) 


Collecting terms and applying trigonometric identities for the difference 


of angics we obtain: 


1 : _ SOStic7 = : ~~ j : S - : 
direction ( ES feeBice? { 1c0s¢ Ge Oe 201) i, 1cos$co (26 6") 
sphere 

= oT aad (2251) 
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Thus frem previously develoned relmtiomshins 2.17 end 2.30, expressions 
for auplitude, phase, and direction of the electric field over thea ref~ 
erence sphere can be conbined to yield the vector field solely as a 
Paictien of the particular reflector under investigation and the choice 
of origin of the reference sphere along the aris. 

While simplification could be performed on the resluting expression, 
the field over the reference sphere is, in itself, only a neans to 
eventually describe the fields over a snaller sphere Surrounding sone 
source region, concentric with the reference sphere, Since the Ssimplifica~ 
ticn would not enhance this procedure, it will not be pursued, Rather, 


Poe resulting expression 


, - Ri nkon Ls edk RC cosé - Ll 
reference a rer int tr tetra tere wasn aK 
sphere Oe hike ORS Val? 


( i_1cos¢sin(26 - 6") — 1,1 cosdcos (26 ~ 6') + ising ) 


Ge 2), 


(where Rw? Row? L. and angle relationships are defined by 2,43, 


Pets 2.33, and 2.40 and are listed for convenience in Section 2,4) 


Jends itself readily to tecniques of machine computation, As a eect. 
the following pages contain graphs representing the field components 
in amplitude and phase for a particular choice of paraneters, The 


are calculated for nornalized anplitude and reflect no variati 
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aszinutheal $6 direction. Thus the components Eo! and Er indicated in the 
graph must be multiplied by cosé ance sind respectively to represent the 
trvué nature of the field ever the entirc sphere. Similarly, the Ea com- 
ponent must be nultiplicd by the variation cosé . The computer program 


developed to calculate the field over the reference sphere is included 


EReappendix II, 
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CHAP Ti 


ELECTROMAGNETIC FIELD EQUATIONS, EXPANSION, AND SOLUTION 


fel in the preceeding chapter expressions for the electric field over the 
reference sphere were determined. In this chapter those fields will be 
matched at the reference sphere to an expansion of spherical harmonics 
following from solution to the vector wave equation Jn spherical coordi- 
nates. These harmonics will be evaluated at a particular radius within 
that of the reference sphere to yield the electric field over the surface 
of that smaller sphere enclosing some source region, These fields will be 
compared to the solution of Ricardi for the same region, Finally, these 


results and the comparison of solutions will be interpreted. 


; 
It may be shown ‘that solution to the scalar Helmholtz equation _ 
Gea Va (3.1) 


for snherical coordinates is given by Separation of variables to be: 


sinmd (3.2) 


m 
= bs (kr) P (cos®) { 


‘ 


hin 


ti 

where bs (kr) is the spherical Besse]. function and p (cose) is the 
essociated Legendre function of the first kind. The choice of cosm¢ 
or sinmd depends upon the variation of the fields at the boundary of 


interest. 


b2 








The fields under investigaticn in chapter 2 were shown to be vector 
fields, and the solution to the vector Helmholtz cquation must be 


considered as: 
(v2 + k2)A @ O (3.3) 


where A is a vector field, representing either the magnetic vector 
Betential or its dual, the electric vector potential. It hes been 

35 7 
shown that a general solution may be formulated by considering a set 
OF fields, both magnetic and electric, which are transverse to the 


radial vector, ie In this case the solutions to 3.3 are found to be: 


=~ 


ye ES eee 
je 
Coe) 
Y= yt lapel 
i 


where Y is the solution to the scalar Helmholtz equation, with the super- 
script indicating that ts solution to the vector Helmholtz equation 


for the magnetic potential A or the @lectric potential F. 


G3 








The most gcneral solution for ficlds in a source-free nemegtnesus region 


4 e ® 3 6 tJ + 
in terms of vector potentials is found to be a linear combination of 


fields with the behavior of equations 3.4, wheres 


E= -y x F = jupA + 1 Ve A) 
Jwe 


(3.5) 


ae; 
" 


Fpmes Gi ieee 1 Ciic 1) 
Ju 


Since the solutions for A and F involve the product of r and ¥, it has 
been found convenient to adopt an alternative definition of the spherical 
Bessel Puneet 


AM (2) 
H (sr) = kr hy Csr) 


2 


A (2) 
The superscript HO indicates we have chosen the spherical Hankel function 


ef the second kind, whase behavior indicates propagation outward from the 
6rigin, as applicable to this problem, Anticipating further, since the 
fields we wish to expand are of the form of equation 2.52, where the ver- 


zation with 4 is given by sing and cos , we choose m = 1 . Thus: 


At?) ie { sin? } 

c - » Y ' 

A. y a H (kr) FEE ee 
(3.6) 

AQ) 4 { sine } 

Be t= : bo H (Cer) P, 6eos8) | oe) 


64 








Expanding, equation 3.5 in component form we obtain: 


a2 

















Peace te es 
r resem (CO mana Jomeats 
jwe oro 
; 2 
Pete eo ee 
Met: —-0@ jwer ora:0 
2 
ee a | nies 
re 68 jJucrsineg ordd 
iG ee as Coal) 
jup dr2 
2 
Hoe eee 
rsind d¢ jwur erdé 
yee 
Hy = Se ee eee 
r 20 qwyrsind grade 


Pweerrins to the cemponent E, in 3.7, we note that $ variation for 


6 


fad 


oF e ¢ e e ¢ 
A and yr is the same. Hencé we chocse the apprepriate variation in 


own 


o> 
equation 3.6. We note further thet the functien P* (cos) is identically 


zero for n = 0. The series need begin for n = ], and may be written: 
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2 ie 


If 
A. = £ a UH_(kr) P_(ces8) cos¢ 
b & nS a n 
n=] 
(3.8) 
= Naan 1 
Be per b. H (kx) Pp (cos®) sind 


Substituting into the appropriate components of 3.7 the electric field 


can be forisulated : 

















n(n+1) ME 
2 ey [ ; 
EY Denn H Csr) P (cos8) cos¢ it 
Seeiee 
i oe 1 
: ng -08¢ 
ony a be H (ce) P (ce 6) cos¢ i, 
rsing 
A(2) j 
1 er : s 
(3.9) tee “2 cos¢ aia Hae e) i, 
jwer dr 6 
A (2) 1 
1 Bee Ole (egaa) ; 
satya bo H (kr) sing n = 
r dé 
A(2) 
1 = en ae dH (kr) 
a. om, sind P_(cos6) no Hy 
jwer sind dr 
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For convenience we define: 





A (2) A (24 e 
2 . net 
soir ae p OOH (ke) 
dr d (kx) 


and since we are concerned with the case of propagation in free space 


Le 
Kia kore 2 
O 0 Oo 


Substituting the relationships from above, multiplying through by j/j = 1, 


and including the summation we form equatien 3.10 from 3,9 : 


e 


Cec) 


ta +a n(ntl) a 1 
E = y} n a3 H Cer) Pr 6cos®) cosé i 


Yr jwer 


At, 1 
BP (ees) ... j. 
a (ae 1 -COS i iz Six cos j 
jta., H Ucr) cosé d_ on mee t sind P feo 6) i 


a 6 sind o_ 


Pole) 1 o 
+ b_H_ (kr) | --cosd oe eee Zi, + @lad d Eeee) 
ren ocean aaa | 0 WS ae 
$1710 di 8 
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Peshback : 
Cc wv OSs ; Tm ‘ 
=: ~~ a _ | cCosne dP (cos8) j.. = singe Weocl a3 
Tid name aa ee aaa ss Se? 6 scaaanantiamammimenean ¢} ; d 
oa do sind 
3.11) 
0 rm ap (cos6). 
= rk mceosnd P_fcosd) ip - sinmd on 7 


cind dé 


Swoscituting 3.11 into 3,10 we can ferm an expression for the tangential 
fields. By exploiting certain orthogonality properties of the tangential 
harmonics, we can derive the constants a and be in 3,10 frem tue fields 
specified at the boundary of interest. In our problem, it will be the 
tangential fields on the reference sphere that will be specified. Note 
Goat trom equation 3.10 the two constants a, and b, are sufficient to 
describe the field vector. It must also be noted that development of all- 
expressions in this chapter has preceeded for the case where the origin 
of the sphere for waich boundary conditions were given coincided with 

the origin of the coordinate system for the spherical harmonics. It now 
becomes necessary to stipulate that the harmonics must be expanded about 
QR Srigin at the center of the reference sphere. To adapt these expression 
we need only to substitute 6' for 0, recognizing that the angle is de~ 
fined as in chapter 2, and for the radius variable r we substitute r', 


e@eercatii.e that it is measured from the prime coordinate origin. 
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From the preceeding we can write fer the tangential field components: 


as C21 ro, 
L gay: ‘ O ; e 
q = = ae ek Gre, ee ae Oa 2 on © > cy) ctgue eel lay ma 
Beang in ee a) De ae r) ote ar nA x) "In 
r n=] 


(3,12 


Coefficients area evaluated by forming the vector dot preduct Ceo ; 
; ] ang 

multiplying each side of 3.12 by sinOd6dd6, integrating over the 

interval O~n and 0-27 on ¢ and 6 respectively, evaluating the fields at 


the reference sphere, radius R', The following properties are listed 


Mest: >” 


a _@ ? Z 
eer 2un (ntl) n=q 
fof sindd6dg = 2dntl 
c° .¢° 
Jq din 
0 n7fq 
o .e = 
oe Cia oe sined6dd 0 forralicn ld 


Applying these properties to equation 3,12 we have: 
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C2mH (kR") UGE Ne 
(3.13) 


on 5-1. ‘ r sd 
(21-11) R Jel on Cha 


SR ns SER ES ee oe 


sing'de'dd 





it 


2mH_ (kR") (n(nt1) 3/2 


where E. is the vector representation of the tangential components cof 


t 
the electric field over the reference sphere at radius R', fo perform 
the integration indicated in equations 3.13 we proceed by expanding 

pae field ie in a fouricr series, We choose a series appropriate for 
the fields over the reference sphere, in this case a scries with even 


symmetry about 6' = 0, and where the 6 dependence matches that ef the 


fields at the boundary, as sugsested by equation 2.52. Procceding: 


(ore, 


2 
wo 


nr 7 ) f ei : en6' , 
Ee CR os 0nd ) fe ( c,, cos Lee d ing ice 08' (3.14) 


The fourier coefficients are obtained in the normal fashion, and are 
covered in appendix 2. 


Consistent with the notation of Ricardi we define coupling coefficients 


1" = f"cosmo!' Paeecorn dé! 
mn O n 
(3.15) 
IL ' 
a = f cosmd' sind! : P coos? ) dQ! 


TOL @) i ; Sos wth 2c ett a 
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Pesctitution of 3.15 into 3.13, using the definition of toarier ce 


Pee icients from 3,14 we obtain: 


o> 


~4 1 
a= 3 es yy (Ce ic =e ae 
n aa Sa m “rn m “mn 
rp? ye 
C3 Gy 
a anoe elite wa 
n seeleeeieaniianient m “nn m nn 


t — 
He (KR ) m=0 


where 


2n + 1 


PAE ID COE 


2(n(n+1))2 





a 


Therefore to fully express the electric field in the spherical harmonic 
expansion we need only te evaluate the fourier coefficients in 3.14 and 
the coupling coefficients given by 3.15. These coefficients, when sub- 
Stituted into the series 3.16 yield ccefficients of the spherical 
harmonic expansion, 
4, 

These constants are normally evaluated by a recursion technique 
rather than solving for each term by evaluation of the integrals. This 
recursion technique is particularly casy to realize with machine computa- 
tion, and is inherent in Ricardi's rethod of evaluation. Thus coefficients 
Will be computed rather than analytically determined. 
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In summary the relationships that express the clectric field as a 
spherical expansion of a spacified field on a given reference sphere 


Such 2s shown in figure 12 are listed: 


ad = EN: A (2) 
HGURENyy) ae es 2y DOH) dex!) Ph(coso') i, 
Gee n= 1. joer. 
(3.17) 
rl Ge) 5 A(z) : e 
Aaa ta . oh c @ & D 
ae Gees 2h f USE ) ae Jt oii pas ) "in 
where oe and ie are given by equation 3.11 with m=1 ; 
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: Sf ad : 
and the coupling coefficients eg Ya piven by 3.15: 
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The fouricr expansion coefficents C. and qd. are defined in equation 3.14: 
fs ; 


oO 

= & sare 9 = so i Sino od sme! 
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Et is the vector representation of the electric field tangential te the 
reference sphere, as determined from 2.52, 

While the series represcntation gives an accurate solution for the 
eu Of all terms, it is recognized that this is impractical to realize. 
Am important section of Ricardi's thesis dealt with justification of an 

; , ea) es 
appropriate truncation of the series solution . That analysis deals with 
Subjects beyond the scope of the present work, although it will be noted 
that the series truncation stems from a desire to avoid generation of 
highly reactive waves that do not, in practice, yield that much accuracy 


to the end result. On the basis of rigorous analysis regacding genera- 


a 
e s 2 4? e e ad 
tion of these highly reactive waves , it will suffice here to state that 
satisfactory results are obtained with the spherical expansion when the 
number of terms of the series is on the order of N = kRy , where Ry is 
the particular radius at which evaluation is desired, 

Therefore, equation 3,17, evaluated at r' = RL , where Ri is the radius 


f f 


enclosing the transmitting feed which will produce the desired field 


E 
reference 
sphere 


Specified by equation 2.52, which will in turn produce after reflection 


a desired field distribution at the aperture, yields 


>? ' 
E eee = ah a we 


sphere 


tims 








The equations for magnetic fields over the feed sphere can be obtained 


from Maxwell's cauations. In particular: 


Hea _ ee Ere al 
sphere jw sphere 


In this analysis we have concentrated on solution cf the electric fields 


at various boundaries. To produce the desired fields at the aperture it 


Mmemiet necessary to separately excite both E and H since 
sia P ad feed i felcue’s 


sphere sphere 
through Maxwell's equations, the existance of one field specifics the 


ether. 


15 








3 


me to. illustrate the principles outlined jp wece cau 


Sections and to show that these principles provide an 
economical and efficient eco of mathematical expansion 
fc telds from the geometrical optTIcs reference sphere to 
an interior feed sphere, a particular case was chosen for 
evaluation, 

Piemearameters JOr batberest Wene wcll occ | Ones ee— mers 
fmeeeesenting a displacement of the origin of the feed and 
f[werence spheres to z = .392R for a reflector of radius 
R = 93 and subtended angle of 37,4 degrees. The field 
misetribution of magnitude and phase were presented in 
eraphical form ine sect tony 225) for vehis chotlee en maweame ecu 
Meese 3 of section 2.1 illustrates the geometry of the _ 
Seeentric coordinate system. 

It must be noted that the choice of eccentric origin 
Orosen for this case of investigation was rather arbitrary, 
femewes the choice of inner sphere radius. This choice was 
Pumemecd, however, from principles of geometrical optics. The 
rays incident upon the extreme edge of the reflector surface 
fweemererlected into the focal region, forming the caustic 
meeace. The points at which the edge rays are tangent to 
mimee caustic surface define the marginal surface focus, that 


mmenme surface generated by the rotation of B-B* of figure 


Mmemobout the axis of revolution of the reflector. 
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Note that all reflected rays must pass through this surface, The 
smallest such surface through which all reflected reys must pass is 
defined by the clircle of least confusion, generated by rotation of 
A-A' of figure 15 about the axis of symmetry. 

While geometrical optics wrongly predicts infinite fields in the 
vicinity of caustics, the fact remains that converging rays indicate 
an area of energy concentration. Hence, an intuitive choice of probable 
ee tion for a minimum-sized feed would be that aperture which most 
closely coincides with the circle of least confusion associated with 
a particular reflector. It is a simple matter to show from the geometry 
that a choice of origin between z = .5R and .6R requires a feed sphere 
eamtne omaer of 5-10 \ for a reflector of 93 X with subtended angle 
37.4 - as considered in this study, to meet these specifications. - 
Hence, a representative distribution of fields over such a feed sphere 
that would be required by the tangential boundary conditions to yield 
a uniform plane wave over the reflector aperture is given by figures 
16, Amplitude and phase of three ee of electuic field*over 
the surface of a feed sphere of radius 8.4 X with origin at z = .52R 
for a reflector of radius R = 93 A are presented = ineiieuccs. |o VA. 
This solution is a graphical result, obtained from evaluation of the 


series solution given by equation 3.16 and 3.17 of section 3.1. 
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Mmepror the choice of parameters indicated an sections. 


the magnitude and phase of the field distribution over the 
memerence aa. as given by equation D5 were presented 
mimeraphical form as figures 13 and 14 of®section 2.5. The 
mearly identical distribution of tangential components and 
the near-vanishing radial component of the field are due 
Bemehe relationship between angular variables 6 and 6 , 


mirere to first order ; 


Only the tangential field components are used in the 
memmmaary condition to generate the fourier coefficients; 
as noted in appendix II. The tangential fields were expanded 


am the following manner from equation 2.52: 


Zz . : ju 
an ( Ecos $i + E SLEO™ 5 el 


Q' > 


3 an = E,cosd (cosh + jsiny) i, + Eee Ceol + join) 2) 
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ihe real ani imaginary parts, dess the cos > sng wdepemdanrcec 
were then plotted graphically and are presented as figures 
meeond 18 . Due to the similarity of magnitudes pesultiange 
Mmeom the nearly-identical nature ye Pemecmeiel components 
[eacated by figures 13 and 14 of section 2.5, only two 
mecher than four plots are presented here. The calculations, 
Meee ver, retain the existence of four separate componenes, 
Pwo real and two imaginary, of tangential fields. These 
memponents were then expanded into a fourier series of the 


feo rin : 


M ee 


= oS ,; 
Es mZo § CO cos¢ igi + d sing i,)cosmé 


where ne Te are now complex fourier cocfifiieitents determined 
fmemoced in appendix II, and M is the limiting value of 


mitemitohest harmonic desired, given, as noted previously, 


by 


<4 
I 

~ 

4K 


Peetl +rithe radius of observation. 
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To check on the accuracy ee the fourler expension technique, the series 
was summed using coefficients determined by the geometrical optics 

meena Over the reference sphere, The results are indicated in fieures 

17 and 18 as small departures from the heavy dark line, which represents 
the valucs used to generate the fourier coefficients. 

These fourier coefficients were then used with the series solution, 
meen by equations 3.16 and 3.17, to evaluate the fields at the reference 
sphere. The trends and accuracy of this series solution, when compared 
to the geometrical optics field initially calculated and subsequently 
used to generate the necessary coefficients, appears to be reasonable 
when viewed in figures 19, It is to be recalled that the series was 
truncated at a value of terms much Jess than infinite in order to avoid 
generation of highly reactive fields, The presentation of figures. 19 
involved 390 terms. As the series’ accuracy improves with the number of 
terms, better agreement between the geometrical optics and resulting 
synthesized fields would theoretically be forthcoming with continued 
summation, It is to be noticéd that the agreement of phase appears 
quite good for the tangential Fields and satisfactory for the radial 
component, Similar agreement occurs in amplitude. The basic boundary 
condition that — erced per nieds saecliits of the tangential fields 
So as to force the field to zero at the reflector surface. Hence, the 
radial component is also determined from this boundary condition, and 
should contain some variation for less than an infinite number of terms. 
The phase of the radial component is similarly determined from the 


real and imaginary parts of the component, and for the radial component 


39 








of near-zero magnitude, the error involved in the difference of small 
numbers begins to override the solution. Hence, for values of radial 
component near zero, rapid phase variations were encountered. These 

were not plotted as not being significant they added no contribution 


to the work. 
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CHAPTER LV 


CONCLUSIONS 


4.1 As noted by the graphical representction of solution offered in 


figures 19, the technique explored by this thesis has yielded results 


ft 


that indicate its accuracy and applicability. A geometrical optics expan=- 
sion from the surface of a spherical reflector to a hypothetical ref- 
erence sphere can be utilizedto generate necessary coefficients required 
by the series solution of spherical harmonics that constitute the 
expression for the electromagnetic fields in the region of concern. 

With regards to the efficiency of the technique, an investigation of 
various eccentric origins from z = .5R to .6R indicated that the radii 
Or the ensuing reference spheres were on the order of Rh’ of, 0R for 
the reflector chosen for investigation by Sag USS Thus, evaluation 
of the series solution at the reference sphere for purposes of compar- 
meeaeco the initially determined geometrical optics field involved a 
reduction in the number of terms that would have been required without 
the eccentric coordinate system by some 40 %. The reduction in terms is 
due to the direct relationship between the radius of observation and the 
number of terms required. Due to the flexible nature of the technique, 
choice of eccentric orlgin to include regions of energy concentration 


near the feed sphere result in the feed radius on the order of 8 A 


for the reflector of 93 \X considered by this study. 





Previous investigation = constrained by a fixed origin and con- 
centric spheres, as given by figure 2, found it necessary to keep the 
radius of the feed sphere on the order of half the radius of the reflec- 
tor. Thus we consider that this technique has provided a substantial 


reauction in the number of terns involved from 


M=kr = ( 2a/d )( 93d /2) 


i 


LO 


Meee) ee 


indicating a savings of a factor of approximately an order of magnitude. 


oe 





mee A next logical step in pursuance of this technique 
would be to concentrate on the series solution at the 
medidt. Of various feed ieee This study suggests that 
meethe fields given by figure 16 can be duplicated over 
M@emsuirtiace of the feed sphere by some “Source within tlre 
moumere the desired fLllumination of uniform plane wave 
meea the aperture of the reflector will be realized. The 
calculation of the fields at the aperture by an alternative 
Memsirement technique, or by experimental methods, if the 
fields over the feéd sphere can be realized by some practi- 
cal feed system, would be a most worthwhile endeavor. 

Further study for a series of parameter variations for 

Sremererttector could, as well, yield a meaningful cee neon 
Beets technique. Based upon the field distribution over 
the feed sphere, design and experimentation with a feed 


system could provide the basis for additional study of this 


mopte at various academic levels of concern. 
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APPENDIX I 


SIMPLIFICATION AND MANIPULATIVE DETAILS OF CERYAIN EXPRESSIONS 


Derivation cf Equation 2.12 


From equation 2.11: 


Ni oes) Ses Vea nil 


Taking the dot product with N to both sides of A.J 


l/p = -NeSx Vx S 


applying the vector identity 

eth x C) = Cath = B) 
then 

Wo = = (ON x Sey x S) AZ 
By cquation 2.10 


S = VL/n_ 


{00 








Substituting and operating: 


VL 1 1 





Vee S = Vox » BN go SVL el re 
Y us ia 
Pee VY x VL = O 0 
a 
V 8 S aa V es x VL 3 aS 
1) 
Pupstitutine A.3 into A,2: 
” @ M 
1/p a ~(N a Sie V ans x Wi) A.4 
) 
We substitute 
v- . *= = sv GE 
is if 
i! 
l/p = -(N x S)e( - V(ln » ) x VL) Aes 
" 
But by the vector identity 
hex Bs ee = CAC) Be) hei c) 


ho] 





it 


a C (SeVin 0) (G-VE) =] Gi-VL Ge ne 
nN 


A.6 1/p 


Fren 2.10, where SeVL 


ae 


y 


H 


1/p = : C NeVin n —- (NeVL)(SeVi1n n) ee: 
n 


The secend term in A.7 is identically zero since by A,.] 
N= = o(S x» V x S) 


ena from 2.10 


then substituting the above into the second tern: 


CNR C85 ox x 5S) 
pr 


Again using the vector identity 


AeCh x C) = Cork x B) 


$02 





Wwe SCC 


Therefore, from A./ 


l/p = NeVin n hoo = 7 le 


Derivation of equation 2,37 ' 


3 


From equation 2,34 


@ x2 ( g2@ = Quz + 22 ) tan“20 fe 9 
a) re) 
From equation 2,35 
R'2 = x2 + 22 ~ 2ecR + c2RA #10 


Supetitution ef A.9 into A.10 
R'2 = ( 22 ~ 227 + ne Jtan220 + z% - 2z2cR + c2R2 
All =(tan226 + 1)z2% ~ 22(z tan? 20 ++ cR)+ z@tan*20-+c*R 
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Multiplying A.ll by cas*Z6 
R'2cos220 = 
' 2 


a ee 
@) 


pelving A.i2 for zg : 


(sin220 - 


sin220 + 


cos*20)22 + 2z(2 sin®20 ‘++ eReos*26) 


c7R2c0s220 A.12 


O- 2 - 22 (2 sin?20 + cReos*20) + Zz Siné 20 ~- c6s226(R'2 = c2R2) A113 


From 2.36: 





a 


3sin26 


Sawetituting A.14 into A.13 


Oe. CRs Oe. re 


sin6d 


sin6 


or Let hase foe te 


R Pie }. dy 


sin2e@ 


g* — 2u(sinOsin20R + cReos220) + sinZ2eR? ~ cos*%20( R'4% = c#R4) = 0 


Ae 15 


1o¥ 





h@emalizing A.15 te the radius ef reflector R 
t 
Z ie 
Zz eee ; 
Th © £202 De yi erscnc ent 
R R 


nA “- az (sinesin2o + c%c0s226) + sin?6 -cos? 


Derivation of equations 26435 2.4% 


Normalizing equation 2.41, 2.42: 


Ry 2 VCR EI ECR ey 
im mn mn 


where 


Zz = cos6 
mn 


x = sind 
mn 


1 { 3ces8 = cos36) 
cn 


j. T, a * 
x # ( 3sin6 = sin36) 


rg 


sind 


On estes ee 


gin20 


Substitution 


05S , 


ve 
cae 


0,A.16=2.37 


fey 





but cos2é@ 


Jin 


2 
in 


then 


in 


SamiLlariy, 


2 
Hoon 


H 


it 


HH 


| 


Qnsagete Oo 


i. ( 2ce@sGeas30 ++ 25inOsin3se + 2 ) 


eva 


ae { 2 + 2cos20 ) 
16 


1 - 2sin26 


4 (1 + sin’6) 
16 





£08? 0 


( cosé ~ y2 + sin6 
sin26 
cos? + sin*9 + Sino” =. —2ebscsrnp 
sin26 51120 
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L Weir! 


Tey ES: 


A,20=2,.41 


oe Z|. 
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COMPUTATIONAL TECHNIQUES RELEVANT TO PROBLEM SOLUTION 


As noted in section 2.5, expressious fer fields weré pregrammed in 
PORTRAN te yield a numerical solutieon to equation 2.52. A cepy ef the 
pregram utilized follows this discussion. Briefly, the legical progression 


Te @Butlined: 


l. For a given reflector of radius R and subtended angle 
Om, equation 2.3] is used to calculate the reference 
sphere radius R' for a particular choice ef parameter 


Ce 


2. Equatien 2.40 is proerammed to yield the angle 6' 
measured from the new origin at z=cR to points on the 
reference sphere corresponding to rays incident upon 


the reflecter surface at a given angle 0, A one-to-one 


correspondence is established between the angles, 


3. Equations 2.41 and 2.42 are programmed along with 


2.33 to yield values for the field anplitule at var- 


\ 


ieus angles 6', measured from the new crigin, on the 
surface ef the reference sphere, The amplitude expression 
is then piven by the inverse distance relationship of 


the form @f equation 2.17. - 


1023 : 





4, The polarization features ccveret by equations 2.46-2,5} 
are programed te indicate the effect 4f redefining the 


ficld cemponents in the new ceerJinate system. 


3, Thé phase cxpreasici is calculatem fer the rcflecter 
undér corsidecratien by relation 2.30 aud 2.33) and, the 
resulting phase is chesen relative te the phase at 


O= O0' = 0 degrees. 


6. The information generated by the previous steps is 
collected and Gutput, if desired, if hade available, 
répresenting & humerical selution te equation 2.52, 


The graphs included in section 2.5 were prepared from results ge 


es) 
ry 
fo 
ct 
ce) 
fs 


as described above, 

it is moted, however, thet the actual field @u the reference sphere 
is only ef passing interest at this point, and that expansion of this 
field into a cesine fourier series, as given by equation 3.14, is ef 


paranount. interest: 


: co 
E reference = 2 (ec cosh i,, + €_ sind i,)cesnd' Aad 3 
ri 9) a > 
sphere n=0 


tangential 
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fhe subroutine FORC was developed to penerata coefficienés fer the 
fourier series from the generated field cenponents. A numerical technique 
fer evaluating the ceoecfficicnts for such a series apprepriave fer 


4d, 
a pericdic sect cf discrete values of a function was available and was 


used as the basis fer this subroutine. As a check, the feurier serics 
Was suitaca, indicating accuracy to within a few percent when viewed 
everall and comparcd te the points used to gencrate the coefficients. 


The technique follows fren consideration of expansion of a periodic 


fumction £(0) in e@ feurier serics; 


£(6) = a /2 + 2, cos760 + b  sinnéd A.24 


where from orthogonality it may casily be shewn that the resultant 


ceef{ficients are given by: 45 
1 21 
= f £(8) cosm6 <6 m=0, 1, cece A. 25 
Ht <= 
i 


As the case of interest involves symmetry sbout 0 = 0, it is obvious 
that a cesine series will suffice to fully express the function and 


Ll will be identically zero fer all un. 
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The integral expression A.25 can be related te a suniatien of discrete 

valucs by appreximations of besic calculus, wnaere the interval fren 

O< 0 < 27 is divided inte 2N intervals, previ@ing ZN + 1 cata pewinte. 
: , =) 


Thus we can formulate a relationship between the number ef intervals 


and the spacing involved where 


2N = 360°/ DEL 


DEL being the spacing between successive points.Thus we approximate, 


fer a sufficiently large number of points: 








é9 = A® = wy ; 
nA a ae | 

6 = nA@ = ot where n= O;1,2, ....2i 
TAGE ieee | 


Substitution of the above approxinatious inte A,25 results in 





2N 
2 MT Yt Z 
Eien y rd Gi _) ¢esn “ue . A.26 
2N + 1 n=0 ye ite oe | 2N + 1 


iif 





Thus we see the functien must be made available in steps @f angle, 


dé@termined by the index n, er at the following discrete points: 


272 
© = i cm OP AL rr eran 


7aT| ae | 


Therefore, the following steps are outlined toe provide an indicatien 
of the progression follevwed through hie subroutine? 

1. Since the serics requires even spacing between the data 
points, it was necessary to interpolate between the points 
generated on the reference sphere by successive steus of 
the angle along the reflector surface. As equation 2.40 
allewed no simple inverse relationship between the two 
angles 6 and 6', small steps of the angle on the reflector 
surface were taken ( spacing = .2 degrees ) to insure the 

-accuracy of linear interpolaticn between the discrete values 


of angle 6' gencrated. 


2. The resulting field expressions were arranged in compenent 
form and subsequeatly expanded frou the exponential represent~ 


ation to the trigonometric forn ;: 


H2 





Le 


(e 
iN 


AC cost ‘+ jsinw ) = AB + 5AT 


The values of field components, twe real and two 
imaginary, were then used to penerate a conplex 
f@uricr coeificicat @f the for 

Cees CR a Gk 


by the technique previously described, 


The answers of interest were these coefficients which, 
when properly inserted into an existing program for 
the series solution to field expressions 3.16 and 3.17. 


the feuricr coefficients c 
- 


a 4 


and Ce are neted in equation 


3.16. 
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